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Ionospheric modification 
Artificial hole injection 
Equatorial Spread F 
.Numerical simulation 


A numerical simulation code for investigating equatorial Spread F in the collisional Rayleigh- 
Taylor regime is utilized to follow the evolution of artificial plasma density depletions injected 
into the bottomside nighttime equatorial F region. The 70 km diameter hole rapidly rises and 
steepens, forming plasma density enhancements at altitudes below the rising hole. The distribution 
of enhancements and depletions is similar to natural equatorial Spread F phenomena, except it 
occurs on a much faster time scale.v 













IONOSPHERIC MODIFICATION: 

AN INITIAL REPORT ON ARTIFICIALLY CREATED EQUATORIAL SPREAD F 


I. Introduction 

Since the discovery of a large dropout (ionospheric hole) in the 
ionospheric total electron content during the launch of NASA's Skylab 
(Mendillo et al ., 1975 a,b) much investigation has gone into the 
controlled artificial production of such holes or plasma density 
depletions in the F region ionosphere (Bernhardt et al .. 1975; 
Bernhardt , 1976; Bernhardt and da Rosa , 1977; Mendillo and Forbes . 
1978; Pongratz et al ., 1978; see also Trans. Am. Geophys. Un ., 58, 
454-455, 1977, Spring Meeting Special Session on "Artificially 
Created Holes in the Ionosphere"). Typically investigators have 
focused on the chemical release of ^0, Hg or high explosives. These 
contain molecular species which react with ionospheric 0 + in such a 
manner that the 0 + loss rate due to these reactions is several orders 
of magnitude faster than natural ionospheric loss rates for 0 + 
(Mendillo et al ., 1975 a,b; Bernhardt , 1976). Recently, Pongratz 
et al . 11978] conducted two high explosive (88kgm mixture of nitro- 
methane and ammonium nitrate) F region chemical releases (260 and 
280 km altitude), codenamed Lagopedo, over the Hawaiian Islands 
(midlatitude) and formed a hole 100 km in diameter) via the 

predicted chemical processes. 

Note: Manuscript submitted March 22, 1978. 
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Within the context of artificial hole formation much of the 
theoretical work to date (Mendillo et al . t 1975 b; Bernhardt , 1976; 
Pongratz et al ., 1978) has centered on the complicated chemical 
reactions and molecular diffusion equations for producing these holes 
(e.g., electrodynamic forces, such as polarization electric fields, 
on the ions or electrons are not considered), lu the past few years, 
experimental ( Kelly et al ., 1976; McClure et al .,1977) and numerical 
simulation studies ( Scannapieco and Ossakow , 1976) have observed 
(under appropriate conditions) the formation of natural, rising 
plasma density depletions (holes or bubbles) in the equatorial 
Spread F ionosphere. In the present paper we wish to perform a 
marriage between the two studies and examine the behavior of a large 
plasma density depletion (hole) artificially injected into the bottom- 
side nighttime equatorial F region ionosphere. 

II. Theoretical Model and Numerical Simulation Studies 

In the present simplified model we will take the approach that 
the large hole has already been created and thus represents a new 
initial condition for our previous nonlinear numerical simulation 
equatorial Spread F studies [ Scannapieco and Ossakow , 1976; Ossakow 
et al ., 1978 a,bj. In this way we do not follow the chemical kinetics 
of the actual hole making process, rather we follow the evolution of 
the hole after it has been produced. Pongratz et al . [1978] noted 
that the entire hole in Lagopedo was produced in less than 5 minutes 
and Sjolander et al . [1978] noted, using rocket in situ ion mass 
spectrometer measurements during Lagopedo, that the dominant ion in 








the artificially created hole was 0 + . Thus, we are able to utilize 
the two dimensional (perpendicular to the magnetic field, B) numerical 
simulation code developed for natural equatorial Spread F studies in 
the collisional Rayleigh-Taylor regime, which is represented by the 
following equations (see Scannapieco and Ossakow , 197b; Ossakow et al ., 
1978 a,b) 

If - f Wi * i) * -V”-" 0 > <» 

7.(v in n7*i) - f (g x z)-Vn < 2 > 

where n is the electron density (n refers to the background equilib¬ 
rium density; n = n Q + n^), is the polorization (induced) potential, 
v R is the usual (used in natural Spread F studies) ionospheric 
recombination rate (see Fig. 1), v in is the usual ion-neutral 
collision frequency, g is gravity (pointing in the -y direction), 
z - B/|B| and V refers to the x and y directions. In our simulations 
realistic altitude profiles of v in and v R (see Fig. 1) and n Q (see 
Fig. 2) have been used. 

The numerical simulation plane coincided with the equatorial 
plane spanning an altitude (y) range from 172 to 452 km in one case 
and 252 to 532 km in the other case (both with Ay - 2 km) with an 
east-west (x) extent of 200 km (Ax ■ 5 km). The system of eqns. (1) 
and (2) was Initialized as follows 







n/n 


exp[o(r-r o -15)^] - 0.97 , 15 km < r-r Q < 35 km 
1.0, r-r Q > 35 km (3) 

where a ■ ln(1.97)/(20)^. This provided a constant 972 depletion 
over a central 30 km of the hole and the depletion decreased to 
ambient density over a 20 km extent on each side (total extent of 
bubble was 70 km). The thickness of the hole at half depletion was 
"v 63 km (see for example Pongratz et al .. 1978). The results of a 
simulation with r Q corresponding to x q m 0, y Q m 300 km, i.e., release 
altitude - 300 km, are shown in Figs. 2-5, which depict contour plots 
of n/n Q (or equivalently n j/ n Q ) at t ■ 0, 500, 1000, and 2000 sec. 

The t ■ 0 frame shows the initial state after the hole has been 
formed by chemical processes. At t * 500 sec the hole has begun to 
steepen on its topside as it rises toward the F peak (altitude * 

354 km). This corresponds to a rise velocity ^ 28 m/sec. We note 
that the initial hole is no longer circular and plasma density 
enhancements are forming to the sides of the bubble (hole). As tne 
initial depletion rises on the bottomside it displaces regions of 
higher density which fall to lower altitudes creating enhancements. 

At t * 1000 sec the top of the bubble is at the F peak and the enhance¬ 
ments are getting larger. This should be contrasted with the natural 
Spread F developing bubble, in this ionospheric profile, which took 
^ 8000 sec to reach the peak (see Scanniipieco and Ossakow , 1976; 

Ossakow et al .. 1978b). At t - 2000 sec the hole is well past the 








peak and has bifurcated on ita topside. Thia is similar to barium 
cloud atriation bifurcation phenomena (aee Zabuaky et al ., 1973; 
Oaaakow et al .. 1977) and would lead to smaller acale structures 
if the present simulation had more resolution. The hole is now 
accompanied by large enhancements and the enhancements in turn are 
forcing depletions to lower altitudes. This picture is far different 
from the initial circular hole at t - 0. Thia figure shows a wide¬ 
spread rapid equatorial Spread F condition. 

Figures b-9 exhibit the results for an artificial depletion 
released at 380 km altitude (x Q ■ 0, y Q - 380 km) in an ionosphere 
with the altitude of the F peak equal to 434 km. In this case the 
time evolution of the hole is much more rapid. For example,between 
t ■ 0 and 50 sec the hole is rising with a velocity •v 88 m/sec. The 
initial rise velocity, here and in the previous case, is consistent 
with a bubble rise velocity given by (Oaaakow and Chaturvedi , 1978) 

Vg “ (g/v^ n ) f (n^/n Q ), where f (n^/n Q ) is an increasing function of 
the percentage depletion, n^/n^, and dependent on bubble shape. In 
particular f corresponds here to a circular-like shape. The case 
depicted in Figs. 6-9 is also to be contrasted to the natural Spread F 
case in this geometry where a bubble took 1000 sec to reach the 
F peak. 

III. Concluding Remarks 

We have presented numerical simulation results, in the 
colllsional Rayleigh-Taylor regime, of artificial depletions injected 
into the bottomside nighttime equatorial F region. It has been shown 
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that such depletions, while they steepen and rise, will cause plasma 
density enhancements to be formed. The entire final picture is 
similar to natural equatorial Spread F numerical simulation phenomena 
except that by artificial injection it occurs on a much faster time 
scale. 
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